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Abstract 



We investigate the properties of large-scale structure predicted in a class of mixed dark 
matter models in which the volatile component (made of particles with high rms velocity) 
derives from the decay of a heavier particle. Such models based on cold+volatile dark matter 
(CVDM) differ from the standard mixture of CDM and massive neutrinos, usually known as 
CHDM, in that they involve a component which has a non-thermal phase space distribution 
function. As a consequence, and differently from CHDM models, the value of the redshift 
at which volatile particles become non relativistic, be varied almost independently 

of the volatile fraction, fix- We compute transfer functions for a selection of such models, 
having 0.1 < fix < 0.5 and different values of z nr . Using linear theory and assuming a scale- 
free primordial spectrum, we compare such models with observational constraints on large- 
scale galaxy clustering and bulk flows, as well as on the abundance of galaxy clusters and 
high-redshift damped Lya systems. We find that these constraints enable us to discriminate 
between different fix and z nr ; within the range of the models inspected, those which can be 
most easily accommodated by the data correspond to the parameter choice fix — 0.2 and 
z nr ~2x 10 4 fl x - 



Subject headings: Cosmology: large-scale structure of the universe - galaxies: clusters: 
general 
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1 Introduction 



The outstanding problem of modern cosmology is to understand the formation of galaxies, 
clusters of galaxies and large-scale structure in the expanding universe. A scenario in which 
these structures form by the growth by gravitational instability of small-amplitude primor- 
dial density fluctuations in a universe dominated by weakly interacting non-baryonic dark 
matter has been the standard framework within this problem has been discussed for many 
years. For most purposes the crucial theoretical quantity in the gravitational instability pic- 
ture is the primordial fluctuation spectrum, P(k). In most theories, density perturbations 
originate as a quantum phenomenon in the very early universe with a power-law spectrum 
of the form P(k) oc k n ; n very close to unity is favoured by many versions of the inflationary 
universe picture, as well as by more general considerations. With the discovery of temper- 
ature anisotropies in the Cosmic Microwave Background Radiation (CMBR) by the COBE 
satellite (Smoot et al. 1992), it is possible to fix the amplitude of this spectrum in a relatively 
unambiguous way on very large scales, of order 1000 Mpc. On smaller scales, the shape of 
the primordial spectrum is expected to evolve from its initial power-law form because of 
the action of various damping and dissipative processes. During several intermediate stages 
different components can have different spectra and, in some models, residual differences can 
still be present at the onset of non-linear stages. 

In many respects the problem of explaining structure formation in the gravitational in- 
stability picture can be reduced to that of finding a power spectrum whose primordial form 
matches the COBE-inferred amplitude on large scales, and whose evolved form simultane- 
ously matches the statistical properties of galaxies and clusters on smaller scales. This has 
proved to be a non-trivial task for the Cold Dark Matter (CDM) model, which, at least in its 
standard formulation (fi = 1, n = 1, h = 0.5 and Gaussian adiabatic fluctuations), is now 
generally accepted to be ruled out by the data (Wright et al. 1992; Taylor & Rowan-Robinson 
1992; Liddle & Lyth 1993; White, Efstathiou & Frenk 1993). The essential problem of this 
model is that, once it is normalized to match the measure CMB temperature anisotropy on 
large scales, it has too large a fluctuation amplitude on small scales ^ 10 h~ 1 Mpc. Despite 
its failure, CDM is nevertheless considered as a reference model, several modifications of it 
having been suggested in order to remedy its shortcomings. Among the "second generation" 
of CDM-based models is the Cold+Hot Dark Matter (CHDM) model, which assumes that 
part of the dark matter content is in the form of massive neutrinos of mass m v ~ 10 eV (Val- 
darnini & Bonometto 1985; Bonometto & Valdarnini 1985; Achilli, Occhionero & Scaramella 
1985; Holtzman 1989; Schaefer, Shafi & Stecker 1992; Schaefer & Shafi 1992; Davis, Summers 
& Schlegel 1992; Holtzman & Primack 1993; Liddle & Lyth 1993; Klypin et al. 1993). In 
this scenario, the small-scale power is suppressed by neutrino free-streaming by an amount 
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which depends on the hot percentage; the number of neutrino species participating in the 
hot component also plays a role (Primack et al. 1995; Babu, Shafer & Shafi 1995). This 
model, with ~ 20-30% to the density from the hot component appears to provide a good 
description of structure on a rather large range of scales, although stringent constraints on 
the exact amount of the hot component are provided by the abundance of high-redshift 
objects (e.g. Ma & Bertschinger 1994; Klypin et al. 1995). 

In this paper we discuss a scenario similar to the CHDM picture, where dark matter has 
a cold component (CDM) and a further volatile component (VDM), which has a phase space 
distribution resulting from the decay of a heavier particle species. A previous paper (Pierpaoli 
& Bonometto 1995, hereafter PB95) calculated the effects of such a particle species upon the 
evolution of fluctuations through the period of recombination and up to the present epoch. 
Results were presented in PB95 in the form of transfer functions for several examples. In the 
following analysis we consider different models with respect to those discussed in PB95, so 
as to provide a wider sampling of the cold+volatile dark matter (CVDM) model parameter 
space. Furthermore, we go beyond the calculation of the transfer functions and submit the 
models to a number of explicit tests, by comparing them with observational data. The 
tests performed take into account: (i) the large-scale galaxy clustering as deduced from 
the analysis of volume-limited galaxy samples obtained from redshift catalogues; (ii) the 
behaviour of bulk velocities calculated using POTENT; (iii) the observed abundance of 
galaxy clusters; (iv) the observed abundance of high-redshift structures traced by damped 
Ly-a systems. Such tests therefore refer to scales ranging from a hundred of Mpc's down to 
a fraction of Mpc. 

The plan of the paper is as follows. In Section 2 we introduce the CVDM models and 
the corresponding linear power spectra. By only resorting to linear-theory approaches, we 
compare in Section 3 such models to observational data. In Section 4 we discuss the main 
results and draw general conclusions from our analysis. 

2 The models 

The model we have in mind is one with two dark matter components, as well as baryonic 
material. To restrict the possible parameter space we consider only models with Q Q = 1, 
Hq = 50 km s -1 Mpc~ x . As for the baryonic contribution, we consider the two values fib = 
0.05 and fib = 0.08, which correspond to the central prediction of standard nucleosynthesis 
and to the 95% upper limit allowed by this constraint (e.g. Reeves 1994) given our chosen 
value of H . 

In conventional Cold+Hot DM models, the hot component is assumed to be made of 
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neutrinos with a mass of a few eV produced in thermal equilibrium in the early universe. 
In such a picture the redshift at which the neutrino becomes non-relativistic is z nr ~ 1.4 x 
10 4 (m !y /10 eV) and the contribution of z/'s to Vl Q is VL U ~ 0.21 g„(m u /10 eV) (g v is the number 
of neutrino spin states with mass m u , originally in thermal equilibrium). Accordingly z nr ~ 
6.7 x lO A VL v /g u and therefore z nr and Q u are not independent, as g v takes only suitable 
discrete values. 

In the scenario we discuss here the VDM particles are produced by the decay of a more 
massive particle and are consequently never in thermal equilibrium. For the purposes of this 
paper, this means effectively that the parameters z nr and Qx, where X denotes the VDM 
species, can be varied independently. 

It is also important to stress that there exist microphysical models where decay processes 
do indeed give rise to a cosmological scenario of the kind we discuss here. One such model 
is described by Bonometto, Gabbiani & Masiero (1994) and discussed also in PB95. In this 
model both supersymmetry (SUSY) and Peccei-Quinn symmetries hold, and it is possible 
that the lightest eigenstate of the ordinary neutral fermion partners in SUSY theories, called 
the neutralino, can decay into the supersymmetric partner of the axion (called the axino). 
The non-thermal axinos produced by the decays would be VDM. Plausible values of the 
model parameters are tlx in the range 0.1 to 1.0 and z nr in the range 10 3 to 10 5 . There does 
remain a constraint between fix and the number of massless neutrinos which is imposed 
by primordial nucleosynthesis considerations (see PB95). We have mentioned this specific 
model merely as an illustration; others can be constructed, but their discussion is beyond 
the scope of this paper. 

Whatever the origin of a non-thermal volatile component, it has a significant effect on 
the nature of the evolved fluctuation spectrum. Our criterion to select models in the £l x -z nr 
parameter space can be sketched as follows. Let g* be the number of helicity states which 
are allowed by nucleosynthesis and g v that associated to neutrinos which are present at the 
nucleosynthesis epoch. Then gx = 9* — 9v is the number of extra spin states to be associated 
to VDM particles. If z nr is the redshift at which VDM becomes non relativistic, then the 
limit for Q x reads 

«*;s^ 9 * (1) 

(see PB95). For each £l x , we choose three different z nr values, namely z nr = 2 x 10 4 f2x, 5 x 
10 4 fi x , and 2 x 10 5 fi X - 

According to Walker et al. (1991), standard stellar light element abundances, combined 
with nucleosynthesis calculations, allow up to g* = 7, a value which is also quite close to 
the 2o" upper bound provided by Olive and Scully (1995). Note, however, that Hata et al. 
(1995) recently gave the more stringent constraint g* < 5 at 95% C.L., which even conflicts 
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with the usual 3 neutrino species. In this paper, we adopt the limit g* = 7. 

The first value of z nr given above corresponds to gx = 5, if eq.([l]) holds as an equality. 
Accordingly, the number of allowed neutrino species is N v = g u /2 = 1. This implies that 
only one massless neutrino is present at nucleosynthesis; the remaining two neutrinos are 
quite heavy and have been already decayed. Taking z nr = 5 x 10 4 Q X instead allows gx = 3 
with N v = 2, while the choice z nr = 2 x 10 5 Q X allows one helicity state for VDM particles 
with N u = 3. 

Models with the first value of z nr would thus be inconsistent with the Hata et al. (1995) 
limits. But note also that, contrary to what one might naively expect, the results we quote 
below for higher values of z nr cannot be used straightforwardly to accomodate this more 
stringent constraint because the shape of the transfer function depends explicitly on the 
number of massless neutrino species. 

As for fix, we choose values in the interval 0.1 < Qx < 0.5, with step 0.1. For each value 
of Qx we allow for both the above baryonic fractions. In Table 1 we list the parameters for 
the resulting 30 models, on which we base the discussion of this paper. Note that, with these 
choices for z nr and fix, differ substantially from those presented in PB95 and probe a more 
carefully chosen part of the parameter space. 

For each model, we follow the evolution of the fluctuations in the baryonic (6b), cold (S c ) 
and volatile (5x) components through equivalence and recombination epochs (the algorithm 
used here has drastically reduced the CPU times with respect to the one used in PB95; 
see however PB95 for its details). The transfer function is, a priori, different for the vari- 
ous components. Since we are not interested here in following the evolution of each single 
component, in what follows we define the overall transfer function as the ratio 

T(k) = 'V/' + /f ^ ■ P) 
(d b + d c + x ) z>Zhm 

Accordingly, T(k) is normalized to unity at scales so large that they enter the horizon 
after recombination. In the above definition, z^or is a redshift at which the smallest scale 
considered is still well above the horizon scale. Numerical values of T(k) are provided by 
the parametric expression 

T(k) = fl + Ec^j , (3) 

where the Cj generally depend on the redshift, due to the residual free-streaming of volatile 
particles. The values of the fitting parameters at z = are given in Table 1, with k measured 
in Mpc -1 . The limiting scale down to which transfer functions are computed is 250 kpc. 
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We will not go into a detailed analysis of the relative behaviour of the models here, 
but it is worth pointing out a couple of important trends. Firstly, at fixed z nr and N v , 
the amount of small scale power relative to large decreases as Qx is increased, due to the 
progressively larger effect of free-streaming. On the other hand, at fixed N v and VL X , the 
power on small scales increases relative to large scales as z nr is increased, since a larger 
z nr corresponds to smaller velocities at the present time, and a consequent reduction of the 
effects of free-streaming. 

Assuming that the COBE signal contains a negligible contribution from gravitational 
waves (see, e.g., Lidsey & Coles 1992), the spectrum normalization can be accomplished by 
matching the observed COBE fluctuation spectrum with its predicted form: 

1 //7n\ 4 r°° 

Ci = (\a lm \ 2 ) = — (-f) J o P(k)jf(kx H )k- 2 dk, (4) 

where ji is a spherical Bessel function and xh = 2c/ H . We take P(k) = Ak (scale-free 
primordial fluctuations: Harrison 1970; Zel'dovich 1972). The spectrum amplitude A is 
estimated by matching the quadrupole (I = 2) value of Qrms-ps = 20/zK (Gorski et al. 1994; 
Wright et al. 1994), with an approximate 10% uncertainty on this value (Gorski, Stompor & 
Banday 1995). The value of A required to normalize each model in this way is also displayed 
in Table 1. 

3 Observational tests 

We concentrate here on comparing analytical predictions of the CVDM models with obser- 
vational constraints on large-scale properties of the density and velocity fields, as well as on 
abundances of cosmic structures. More refined investigations involving the use of N-body 
simulations to account properly for effects of non-linear clustering, are beyond the scope 
of this paper. All the results in the following depend only on the linear power-spectrum. 
Results for the CVDM models are also compared to the predictions of the standard CDM 
and the CHDM model with VL U = 0.3 for the fractional density contributed by one massive 
neutrino species. 

3.1 Power— Spectra 

In order to compare the power spectra of the CVDM models to that of the galaxy distribution, 
we plot in Figure 1 P(k) for models with Qx = 0.1, 0.2, 0.3 and 0.5, as well as for two volume- 
limited subsamples of the CfA2 and SSRS2 redshift surveys (da Costa et al. 1994). The 
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spectra plotted are all normalized to COBE, so that galaxies are assumed to trace exactly the 
DM distribution (allowing for linear biasing would just result in a vertical shift of the curves). 
Linear-theory redshift space distortions are introduced by multiplying P(k) by the correction 
factor / = 1 + 2/(36) + 1/ (5b 2 ) (Kaiser 1987), where we assume b = 1 for the biasing factor. 
Since we plot linear spectra, we should take in mind that the comparison with data is reliable 
only for k^, 0.2 (/i _1 Mpc) _1 . At smaller scales, non-linear gravitational clustering makes 
both the linear spectrum shape and the treatment for redshift -space distortions inadequate. 
In any case, we expect the net result of non-linear effects to be an increase of P(k). Those 
models whose linear P(k) already fall above the data points would therefore have even harder 
time if a proper non-linear treatment were performed; this is the case for CDM and for all 
the Qx = 0.1 CVDM models. As Qx is increased, we note that the models with large 
values of z nr display more power than the data for 0.1^ k^, 0.2 (/i _1 Mpc) _1 . This is just 
the consequence of the fact that, at larger values of z nr the volatile particles have smaller 
velocities and, therefore, fall into the CDM potential wells at an earlier epoch. 

In order to provide a more quantitative description of the clustering, we compute a$, 
which is defined as the rms fluctuation amplitude, 



within a top-hat window, W(kR) = 3(sin kR — kR cos kR)/(kR) 3 , of radius R = 8/i Mpc. 
The <7g values for all the models are listed in Table 2 and are compared with that, a§ = 
0.90 ± 0.05, reported by Loveday et al. (1995) for Stromlo-APM galaxies in real space. 
It turns out that models with Q x = 0.1 or with large z nr have the unpleasant feature of 
implying a substantially antibiased galaxy distribution, in some cases comparable to that of 
CDM. A constraint on a% for the DM distribution comes from the cluster abundance, which 
suggests og ~ 0.6 for Qq = 1 models (e.g. White, Efstathiou & Frenk 1993; see below for a 
more detailed discussion of the cluster abundance predictions for our models). We note that 
only small variations of erg are obtained by increasing the baryonic contribution from 5% to 



Smaller values of <j$ may be allowed by the COBE normalisation if one tilts the primordial 
spectral index ton < 1, also possibly allowing for the presence of some tensor contribution in 
the CMB temperature anisotropies. For instance, taking n = 0.9 for the model 4 (Qx = 0.2, 
z nr = Ax 10 3 ) one gets a 8 = 0.70. 

As a further characterization of the power-spectrum shape, we computed the parameter 
r defined by V = 0.5(3.4<7 25 /0.95<7 8 ) ~ 1/a3 (Wright et al. 1992; Efstathiou, Bond & White 
1992); T ~ 0.25 is required by the data. According to the values reported in Table 2, we 
confirm the visual impression obtained from Figure 1: models with Qx = 0.1 have a too large 




(5) 



8%. 
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r (they are too similar to CDM), while models with a larger amount of volatile component 
fare better, unless one takes a large value of z nr . 



3.2 Bulk Velocities 

The rms bulk velocity, Vb u ik(-R)> is defined as the rms matter velocity after smoothing over 
a volume of size R. For Q a — 1 it is connected to the power-spectrum according to 



where W(kR) is the window function specifying the shape of the smoothing volume. By 
comparing eqs.([5|) and (|]), it is clear that Vb U ik(-R) gives more weight to long wavelength 
modes than than or. Therefore, we expect bulk velocities on large scales to depend only on 
the Qrms-ps normalization and not on the profile of the transfer function. 

Reliable Hulk data for top-hat spheres centered on the Local Group on scales of few 
tens of Mpcs are provided by the POTENT reconstruction method (e.g. Bertschinger et 
al. 1990; see also Dekel 1994, and references therein). In Figure 2 we compare our model 
predictions with the latest POTENT data (courtesy of A. Dekel). In order to account for 
the velocity smoothing procedure in the reconstruction method, we convolved the power- 
spectrum in eq.(0) with a Gaussian filter of radius Rf = 12 /i _1 Mpc. As expected, any 
difference between models at large scales is negligible and all of them are in remarkable 
agreement with data. On smaller scales (^ 40 /i _1 Mpc) the predicted Vbuik values tend to 
be larger than the observational one. In this respect, models with Qx > 0.4 perform better, 
thanks to the steep P(k) profile at large k, although this is unfortunately inconsistent with 
galaxy clustering data. On the other hand, independent estimates of bulk flows (e.g. da Costa 
et al. 1995) agree with the POTENT one only on scales R ~ 50-60 h~ 1 Mpc. Therefore, we 
do not regard this marginal discrepancy as a serious problem for any of the models we have 
considered. 

3.3 Cluster abundance 

According to the standard Press & Schechter (1974) approach, the number density of col- 
lapsed structures arising from Gaussian initial fluctuations and having mass larger than M 




(6) 



is given by 




(7) 
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Here n(M) dM is the number density of objects with mass in the range [M, M + dM] and is 
related to the power-spectrum according to 



n(M) dM 



a/27T / JR Or 




— exp 
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where 
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k 4 P(k) 



dW 2 (kR) dk 
d(kR) kR ' 



(9) 



In the above expressions, we assume that the mass scale M is related to the length scale R 
according to M = fpR 3 , with / the "form factor", which is specified by the shape of the 
filter W and p the average density For the Gaussian window, that we assume here, it is 
/ = (2-7r) 3 / 2 , while f = An/3 for a top-hat window. The parameter 5 C is the critical density 
contrast, which represents the threshold value for a fluctuation to turn into an observable 
object, if evolved to the present time by linear theory For a top-hat spherical collapse one 
has 5 C = 1.68, but the inclusion of non-linear effects, as well as aspherical collapse, may 
lead to a lower value of S c . For example, Klypin & Rhee (1994; KR94 hereafter) found that 
the cluster mass function in their CHDM N-body simulations is well fit by eq.(^) by taking 
8 C = 1.5 for the Gaussian window. In order to account for the rather poor knowledge of S c , 
we prefer to compute iV(> M) for different values of this parameter in the range [1.4, 1.7]. 

The results of this analysis are reported in Figure 3, where we also compare them with 
observational results. Values of N(> M) for 5 C = 1.5 are also listed in Table 2. Following 
White et al. (1993), we take M = 4.2 x 10 14 /i _1 M Q for the limiting mass at which to estimate 
the mass function. We prefer not to consider a larger value, ~ 10 15 h~ 1 M & , since this would 
correspond to the exponential tail of the cluster mass function (e.g. Bahcall & Cen 1992) 
and, as a consequence, large variations in the cluster abundance would be associated with 
uncertainties in the cluster mass estimates. The dashed band in Figure 3 corresponds to the 
range between the observational result of White et al. (1993; lower limit) based on X-ray 
data, and that of Biviano et al. (1993; upper limit) based on velocity dispersions. We note 
that realistic observational uncertainties are probably larger than the difference between such 
two results. They may be due to systematic effects, related to assumptions used to connect 
X-ray temperature and DM potential profiles, or to biases in estimating cluster masses from 
internal velocities under the virial assumption. 

Even bearing such warnings in mind, it seems difficult to reconcile with the data those 
models which overproduce clusters by one order of magnitude or more, for any value of 5 C 
value. This is the case for Qx = 0.1 and, in general, for those models having a large value 
of z nr . Even though taking the larger baryonic fraction decreases the cluster abundance, 
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its effect is nevertheless only marginal for those models which have an exceedingly large 
N(> M). 

It is not clear whether such large discrepancies may be overcome on the ground of ob- 
servational biases. For instance, let us consider the model 5 (Qx = 0.2, z nr = 10 4 ) as a case 
providing a large N(> M). If we allow for an underestimate of cluster masses by a factor 2 
(i.e., M = 8.4 x 10 14 /i -1 M Q ; see, however, Evrard, Metzler & Navarro 1995, for arguments 
in favour of precise mass determinations from X-ray data) and take 5 C = 1.5, it would give 
N(> M) ~ 1.2 x 10~ 5 (/i~ 1 Mpc)~ 3 , which is still quite far from the observational band. 

On the other hand, we do not believe that the observational situation is clear enough to 
rule out at a large confidence level models which are discrepant by a factor 2-3 with respect 
to the reported abundances. Even adopting such a rather conservative position, it is fair to 
say that only models with Qx > 0.2 and low z nr are not excluded by this constraint. In 
this respect, the availability of more and more determinations of cluster masses based on 
the independent technique of weak gravitational lensing (e.g. Squires et al. 1995) will be 
extremely welcome. 

3.4 High— redshift behaviour 

A further constraint on power spectra comes from observations of high-redshift objects. The 
most reliable such constraint concerns the abundance of damped Ly-a systems (DLAS). 
These are observed as wide absorbtion troughs in quasar spectra, due to a high HI column 
density (> 10 20 cm -2 ). The fact that at z ~ 3 the fractional density of HI gas associated 
with DLAS is comparable to that contributed by visible matter in nearby galaxies, suggests 
that DLAS actually trace a population of collapsed protogalactic objects (see Wolfe 1993, 
for a comprehensive review). Lanzetta et al. (1995) and Wolfe et al. (1995) presented data 
on DLAS up to the redshift z ~ 3.5, while the recent compilation by Storrie-Lombardi et 
al. (1995) pushed this limit to z ~ 4.25. Based on these data, the latter authors claimed 
the first detection of a turnover in the fractional density, Q g , of neutral gas belonging to the 
absorbing systems at high redshift. 

Several authors recognized DLAS as a powerful test for DM models, based on a linear 
theory approach, Subramanian & Padmanabhan (1994), Kauffman & Chariot (1994), Mo & 
Miralda-Escude (1994) and Ma & Bertschinger (1994) concluded that the standard CHDM 
scenario with Q u = 0.3 is not able to generate enough collapsed structures at zk> 3, due to 
the lack of power on galactic scales. However, either lowering Q u to about 0.2 (Klypin et 
al. 1995) or 'blueing' the primordial spectrum, Pi(k) oc k n to n ~ 1.2 (Borgani et al. 1995) 
keeps CHDM into a better agreement with data. Katz et al. (1995) resorted to numerical 
simulations of DLAS and found that even the CDM model with a normalization as low as 
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o"8 = 0.7 satisfies the DLAS constraint. 

In order to connect model predictions from linear theory with observations, let 

n »" (M -^ erfc (7ffe)' (10) 

be defined as the fractional density contributed at the redshift z by collapsed structures of 
mass larger than M. Accordingly it is Q g = a g QbQ co \\, where a g is the fraction of HI gas 
which is involved in the absorbers. 

One expects the value of a g , to decrease well below unity at low redshift, due to gas 
consumption into stars. Note that varying this number turns into a linear rescaling of 
Q g . Since we assume here that a g — 1, we compare data and CVDM predictions at the 
highest redshift allowed by the data, z = 4.25 (Storrie-Lombardi et al. 1995). We estimate 
<t(M, z = 4.25) for a Gaussian window, by explicitly computing the transfer function at 
that redshift, so as to take into account effects of residual free-streaming of the volatile 
component. 

As for the value of the critical density contrast S c , results based on N-body approaches 
by Efstathiou & Rees (1988) and by Klypin et al. (1995) suggests that 1.3^ 5 C ^ 1.5 for 
a Gaussian window, while Ma & Bertschinger (1994) found indications for S c ~ 1.7-1.8 
for a top-hat window. In the following we will report results for 5 C = 1.5 and gaussian 
window; while the effect of varying this parameter is discussed in more detail by Borgani et 
al. (1995). Lacey & Cole (1994) have realized a detailed test of the Press-Schechter theory 
against scale-free N-body simulations, by checking the effects of varying the window and the 
halo identification method. As a result, they found that 5 C ~ 1.3 and 5 C ~ 1.8 are in general 
adequate to describe the halo mass function for Gaussian and top-hat window, respectively. 
In any case, we verified that rather similar results are obtained either using 5 C = 1.5 with a 
Gaussian window, or 5 C — 1.7 with a top-hat window. 

The results of our analysis are reported in Figure 4, where we plot the neutral gas 
fraction associated to DLAS, for all the models with £l b = 0.05 and compare them with 
the observational data. In the light of all the above-mentioned uncertainties in realising 
such a comparison, we prefer here to adopt a conservative approach and to consider in this 
comparison the result of Storrie-Lombardi et al. (1995; Q g = 2.2 ±0.5 at z = 4.25 for Q = 1 
and h = 0.5) as a lower bound. Consistently, the dashed areas in Figure 4 are delimitated 
by the above la lower limit. Only model falling below this limit are ruled out. The effect of 
varying the limiting mass of the protogalaxy hosting DLAS by an order of magnitude may 
be judged by comparing open and filled dots, which correspond to M = 5 x 10 9 /i _1 M Q and 
M = 5 x lO lo /i~ 1 M , respectively. Numerical values of Q g for all the models are reported in 
column 5 of Table 2, where also results for CDM and CHDM are given. 
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As expected, taking lower fix and larger z nr make easier the agreement with data. All 
the models with fix < 0.2 are able to pass the DLAS test, while larger volatile fractional 
densities are only allowed for large z nr , which however turns into a wrong power-spectrum 
shape. 

Although several models can clearly be ruled out already at this level on the ground of 
DLAS data, nevertheless it is clear that more precise conclusions can only be reached with a 
better knowledge of the variables entering in the Press-Schechter prediction for fi 9 (i.e. the 
mass M, and the parameters S c and a g ). A more accurate definition of what is a DLAS in a 
given DM model can only be achieved with numerical simulations involving hydrodynamics 
(Katz et al. 1995), which would be able to trace the history of galaxy formation. As for 
observations, the possibility that systematic biases may affect the final results have been 
recently suggested. For instance, Bartelmann & Loeb (1995) have recently pointed out 
that amplification biases due to gravitational lensing of QSOs by DLAS could led to an 
overestimate of Q g , by an amount which however decreases with redshift. Fall & Pei (1995) 
argued that dust obscuration may act in the opposite direction so as to bias downwards 
the estimated fi 9 . Verifying the actual relevance of such effects surely requires a substantial 
investment of observational and theoretical effort. 

4 Discussion 

The results we have presented demonstrate that the CVDM hypothesis yields potentially 
interesting models of structure formation. The aim of this paper is to show that rather 
slight changes in the parameters of volatile dark matter can make a significant difference to 
the transferred power spectrum. This contrasts with the case of a cold component, where the 
physical properties of the candidate particle do not really matter at all, in that the physical 
origin of the hot particles may leave a detectable imprint in the clustering pattern. In this 
context it is important to verify up to which point the shape of the distribution function 
causes differences compared to the standard scenario based on relic thermal neutrinos. 

The ability to change z nr almost independently of fix is especially significant in this 
respect: the power spectra we have obtained display considerable variations at a fixed value 
of fix- Indeed, although the CVDM class of models involves one more parameter than is the 
case for CHDM, we have shown that observational data nevertheless allow us to put rather 
stringent constraints on the permitted values of z nr and fix, even at the level of linear- 
theory. The most stringent of these constraints comes from the simultaneous requirement 
for a model to satisfy the observed abundance of high-redshift DLAS and of galaxy clusters. 
As for DLAS, the rather large value of the HI gas fraction involved in the absorbing systems, 
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tig, implies a substantial amount of power on galaxy scales, so as to favour models with 
Qx~ 0.2. A larger volatile component would be allowed only resorting to a high value of 
z nr ~ 2 x 10 5 fix (cf- Figure 4). On the other hand, models with small Vt x and/or large 
z nr behave too much like the standard CDM model, drastically overproducing clusters (cf. 
Figure 3). 

Therefore, the overall result would be that models with z nr ^ 5 x 10 4 f2x have a hard 
time, quite independently of Qx- Among the models inspected, the only model which passes 
all the tests, or at least which can not be confidently ruled out, is the one with £l x = 0.2 
and z nr = 4 x 10 3 . It is worth recalling, however, that such a model with low z nr requires 
that volatile particles occupy at least 5 helicity states [cf. eq.(l)]. We recall that this can 
accommodated only if (a) g* = 7 is allowed by standard nucleosynthesis and (b) two neutrino 
species are sufficiently massive that they have already decayed at the nucleosynthesis epoch. 

An alternative possibility, holding if the physics of the decay is quite different from the 
axino model suggested in PB95, is that the decay itself takes place after the nucleosyn- 
thesis epoch. This would make low z nr . models compatible with all N v . It must, however, 
be remembered that changing N v itself causes an alteration of the transfer function, and 
a straightforward extrapolation of the above results to greater N v values is not allowed. 
Changing the relation between N u and z nr opens the way to inspecting different models and, 
in this context, we should also bear in mind that our analysis has been based on assuming 
a scale-free primordial spectrum, while variations around this model are allowed by some 
classes of inflationary schemes. For instance, taking Pi(k) oc k n with n < 1 (Adams et al. 
1993; Liddle & Lyth 1993 and references therein) decreases the amount of power on the clus- 
ter mass scale, so as to alleviate the problem of cluster overproduction displayed by "colder" 
models. However, the amount of this tilt can not be too large, in order not to conflict with 
CMB (Bennet et al. 1994) and large-scale peculiar motions constraints (Tormen et al. 1993). 
On the other hand, the case of "antitilting" , with n ~ 1.2, has been recently advocated to 
alleviate some of the problems of the CHDM scenario (Dvali, Shafi & Schaefer 1994; Lucchin 
et al. 1995). However, the subsequent increase of power on small scales goes in the unde- 
sired direction as far as the cluster abundance is concerned (Pogosyan & Starobinsky 1995; 
Borgani et al. 1995). 

As a final remark, we should stress that the analysis presented in this paper is only 
preliminary and is entirely based on linear calculations. In order to be more definitive we 
would like to extend it in two main directions. Firstly to calculate more detailed properties of 
the CMBR fluctuations they produce: we anticipate a rather different signature on angular 
scales around a degree than in the standard CHDM models. Furthermore, we would also 
like to study the non-linear evolution of some of these model by performing numerical 
calculations using N-body and other procedures. 
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Figure Captions 



Figure 1. Comparison between observational and CVDM linear power-spectra in redshift 
space. Open and filled triangles are the power-spectrum for two volume limited subsamples 
of the combined CfA2+SSRS2 survey (da Costa et al. 1994). Each panel refers to a fixed 
value of the volatile fraction fix, while the dotted, short-dashed and long-dashed curves 
correspond to different z nr . For reference, we also plot the CDM (solid curves). All the 
models are for fl b = 0.05. 

Figure 2. Comparison between the POTENT bulk flow (filled dots) and that provided by 
the CVDM models. Different panels are for different fix values. In each panel, different 
curves refer to the different values for z nr . 

Figure 3. The abundance of galaxy clusters with M > 4.2 x lO 14 /i _1 M . The shaded 
area is delimited by the observational results by Biviano et al. (1993; upper limit) and by 
White et al. (1993; lower limit). In each panel, corresponding to different fix values, dotted, 
short-dashed and long-dashed curves are for the three different values of z nr . Each pair of 
curves correspond to the two values of fib, the lower one being for fib = 0.08 and the higher 
one for fl b = 0.05. 

Figure 4. The fractional density of neutral gas involved in DLAS at redshift z = 4.25. The 
shaded area is the observational constraint and is delimited from below by the la lower limit 
by Storrie-Lombardi et al. (1995). Each panel refers to a fixed Qx value and reports the 
predicted fl g as a function of z nr . Filled and open dots correspond to M = 5 x 1O 9 /i _1 M 
and M = 5 x 10 10 /i~ 1 M o for the limiting mass of the protostructures hosting DLAS. We 
assume a Gaussian window with S c = 1.5 and a g — 1 for the HI gas fraction involved in the 
absorber. 
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Table 1: Model parameters and power spectra. Column 2: number of massless neutrino 
species; Column 3: volatile fractional density; Column 3: redshift at which the volatile 
component becomes non-relativistic (in units of 10 4 ); Columns 5 to 9: fitting parameters of 
the transfer functions [see eq.flU)]. 



Model N u fix z nr /10 4 A(10 7 ) Ci c 2 c 3 c 4 



n b = 0.05 



1 


1 


0.1 


0.2 


1.300 


-0.6133E+00 


0.1494E+02 


0.1220E+03 


0.6733E+02 


2 


1 


0.1 


0.5 


1.306 


-0.4284E+00 


0.1212E+02 


0.9742E+02 


0.6607E+02 


3 


3 


0.1 


2.0 


1.311 


-0.3665E+00 


0.1542E+02 


0.7916E+02 


0.8917E+02 


4 


1 


0.2 


0.4 


1.318 


0.5110E-01 


-0.1016E+01 


0.1924E+03 


0.1559E+03 


5 


1 


0.2 


1.0 


1.319 


0.1032E-01 


0.4255E+01 


0.1113E+03 


0.1629E+03 


6 


3 


0.2 


4.0 


1.327 


0.4886E-02 


0.1624E+02 


0.3607E+02 


0.1997E+03 


7 


1 


0.3 


0.6 


1.354 


0.1263E+01 


-0.2246E+02 


0.2301E+03 


0.3903E+03 


8 


1 


0.3 


1.5 


1.303 


-0.5240E+00 


0.1509E+02 


0.2086E+02 


0.4061E+03 


9 


3 


0.3 


6.0 


1.266 


-0.1997E+01 


0.5180E+02 


-0.1305E+03 


0.4194E+03 


10 


1 


0.4 


0.8 


1.325 


0.1052E+00 


0.8104E+01 


-0.1939E+02 


0.1053E+04 


11 


1 


0.4 


2.0 


1.248 


-0.2475E+01 


0.5754E+02 


-0.2466E+03 


0.9205E+03 


12 


3 


0.4 


8.0 


1.168 


-0.5488E+01 


0.1240E+03 


-0.4753E+03 


0.8259E+03 


13 


1 


0.5 


1.0 


1.249 


-0.3115E+01 


0.1033E+03 


-0.7761E+03 


0.2781E+04 


14 


1 


0.5 


2.5 


1.107 


-0.7916E+01 


0.1846E+03 


-0.9825E+03 


0.2091E+04 


15 


3 


0.5 


10.0 


1.063 


-0.9233E+01 


0.1943E+03 


-0.8273E+03 


0.1325E+04 



n b = 0.08 



16 


1 


0.1 


0.2 


1.323 


-0.1009E-01 


0.8960E+01 


0.1418E+03 


0.6768E+02 


17 


1 


0.1 


0.5 


1.319 


-0.2039E-01 


0.5736E+01 


0.1202E+03 


0.6457E+02 


18 


3 


0.1 


2.0 


1.318 


-0.6469E-01 


0.6622E+01 


0.1021E+03 


0.8145E+02 


19 


1 


0.2 


0.4 


1.308 


-0.4181E-01 


-0.9723E+01 


0.2449E+03 


0.1497E+03 


20 


1 


0.2 


1.0 


1.348 


0.9112E+00 


-0.1115E+02 


0.1681E+03 


0.1530E+03 


21 


3 


0.2 


4.0 


1.317 


-0.1257E+00 


0.1017E+02 


0.5465E+02 


0.1926E+03 


22 


1 


0.3 


0.6 


1.443 


0.3845E+01 


-0.5984E+02 


0.3617E+03 


0.3668E+03 


23 


1 


0.3 


1.5 


1.346 


0.7783E+00 


-0.3884E+01 


0.8075E+02 


0.3960E+03 


24 


3 


0.3 


6.0 


1.203 


-0.3672E+01 


0.6458E+02 


-0.1643E+03 


0.4436E+03 


25 


1 


0.4 


0.8 


1.392 


0.2294E+01 


-0.3311E+02 


0.1601E+03 


0.1005E+04 


26 


1 


0.4 


2.0 


1.244 


-0.2619E+01 


0.6079E+02 


-0.2866E+03 


0.1030E+04 


27 


3 


0.4 


8.0 


1.159 


-0.5641E+01 


0.1198E+03 


-0.4603E+03 


0.8270E+03 


28 


1 


0.5 


1.0 


1.279 


-0.2064E+01 


0.8336E+02 


-0.6956E+03 


0.2808E+04 


29 


1 


0.5 


2.5 


1.148 


-0.6521E+01 


0.1636E+03 


-0.1016E+04 


0.2456E+04 


30 


3 


0.5 


10.0 


0.925 


-0.1341E+02 


0.2297E+03 


-0.9345E+03 


0.1439E+04 



Table 2: Statistical properties of large-scale structure. Column 2: r.m.s. fluctuations within 
a top-hat sphere of 8 ft, _1 Mpc radius. The observational result refers to the APM galaxy- 
distribution in real space (Loveday et al. 1995). Column 3: 'extra power' parameter (see 
text); observational result from Peacock & Dodds (1994). Column 4: number density of 
clusters with mass larger than M = 4.2 x 10 14 /i~ 1 M o (in units of 1(T 6 ( /i^Mpc)" 3 ). Lower 
and upper values for the observational result are from White, Efstathiou & Frenk (1993) 
and Biviano et al. (1993), respectively. Column 5: fractional density of neutral gas within 
collapsed structures of mass 5 10 w h~ 1 M Q at redshift z = 4.25, in units of 1CT 3 ; the Gaussian 
window with 5 C = 1.5 is assumed; observational result from Storrie-Lombardi et al. (1995). 



Model a 8 T N(> M) Q g 

Observ. 0.90 ± 0.05 0.25 ± 0.05 (4^6) 2.2 ± 0.5 



Vl b = 0.05 



1 


1.05 


0.30 


36 


8.7 


2 


1.23 


0.34 


53 


11 


3 


1.23 


0.35 


53 


10 


4 


0.87 


0.19 


19 


1.4 


5 


1.10 


0.25 


40 


3.2 


6 


1.23 


0.33 


53 


5.1 


7 


0.76 


0.13 


11 


4E-02 


8 


1.02 


0.21 


33 


0.6 


9 


1.22 


0.36 


52 


2.3 


10 


0.68 


0.10 


6.1 


4E-05 


11 


0.96 


0.18 


28 


2E-03 


12 


1.29 


0.53 


60 


0.8 


13 


0.62 


0.08 


3.4 


IE- 11 


14 


0.95 


0.21 


27 


1E-04 


15 


1.27 


0.59 


58 


0.3 


CDM 


1.33 


0.47 


62 


30 


CHDM 


0.86 


0.16 


18 


1E-02 
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Table 2: continued 



Model 


08 


r 


iV(> M) 




Observ. 


0.90 ±0.05 


0.25 ±0.05 


(4-6) 


2.2 ±0.5 






tt b = 0.08 






16 


1.03 


0.28 


34 


6.0 


17 


1.20 


0.31 


50 


8.8 


18 


1.26 


0.33 


55 


8.9 


19 


0.81 


0.16 


14 


0.4 


20 


1.05 


0.21 


36 


1.9 


21 


1.26 


0.30 


55 


4.3 


22 


0.74 


0.10 


9.1 


3E-03 


23 


1.01 


0.18 


32 


0.2 


24 


1.21 


0.35 


51 


1.8 


25 


0.64 


0.08 


3.9 


2E-07 


26 


0.95 


0.17 


27 


4E-03 


27 


1.27 


0.48 


59 


0.8 


28 


0.60 


0.07 


2.5 


2E-15 


29 


0.93 


0.14 


24 


2E-06 


30 


1.19 


0.53 


51 


0.1 


CDM 


1.28 


0.44 


58 


20 


CHDM 


0.82 


0.15 


15 


6E-03 



20 
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